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Abstract We discuss the energy dependence of hadron production in relativistic nucleus- 
nucleus collisions. Several 'anomalies' in the energy dependence have been 
predicted as signals of the deconfinement phase transition and three of these 
signals are observed at the CERN SPS indicating that the onset of deconfinement 
in Pb+Pb collisions is located at about 30 A-GeV. 
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1. Introduction 

The data on heavy nucleus-nucleus (A+A) collisions suggested that there is a 
significant change in the energy dependence of the pion and strangeness yields 
which is located between the top AGS (1 1 A-GeV) and SPS (158 A-GeV) en- 
ergies. Based on the statistical approach it was speculated that the change is 
related to the onset of deconfinement at the early stage of A+A collisions, and a 
simplified quantitative model was developed, the Statistical Model of the Early 
Stage (SMES) [1]. The SMES predicted a sharp maximum in the multiplicity 
ratio of strange hadrons to pions at the begining of the transition region. This 
prediction triggered a new experimental programme at the SPS - the energy 
scan programme - in which the NA49 experiment recorded central Pb+Pb col- 
lisions at several energies: the results from the run at 40, 80 and 158 A-GeV are 
published [2], the results from the 30 A-GeV run were shown for the first time 
at the conference SQM03 [3], the data at 20 A-GeV are still being analysed. 
The energy scan program at the CERN SPS results an observation of several 
'anomalies' in the energy dependence of hadron production in the same do- 
main of the collision energy. These 'anomalies' are interpreted as signals of 
the deconfinement phase transition in A+A collisions. In this report we review 
the physical arguments which lead us to the proposed signals as well as discuss 
their experimental status. 
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2. Signals of Deconfinement: Model Predictions 

Originally, two signals of the deconfinement phase transition were proposed 
within SMES [1]: the energy dependence of the mean pion and strangeness 
multiplicities. Recently, two new signals were suggested within SMES: the 
energy dependence of the inverse slope of the transverse mass spectrum of 
kaons [4] and the energy dependence of properly filtered multiplicity fluctua- 
tions [5, 6]. 

An exact nature of the deconfinement phase transition is still debated. On 
the other hand, it is rather well established in the lattice QCD at zero bary- 
onic chemical potential that very strong changes of the energy density e take 
place in a narrow temperature interval AT = 5 10 MeV. Within this tem- 
perature interval the energy density changes by about an order of magnitude, 
whereas the pressure remains approximately unchanged. One may refer to this 
temperature interval as a 'generalised mixed phase'. 

The Pion Multiplicity. The majority of all particles produced in high energy 
interactions are pions. Thus, pions carry basic information on the entropy cre- 
ated in the collisions. On the other hand, the entropy production should depend 
on the form of matter present at the early stage of collisions. A deconfined mat- 
ter is expected to lead to the final state with higher entropy than that created by 
a confined matter. Consequently, it is natural to expect that the onset of cre- 
ation of a deconfined matter should be signalled by an enhancement of the pion 
production. Clearly, a trivial dependence of the pion multiplicity on the size of 
colliding nuclei should be removed and thus a relevant observable is the ratio 
of the mean pion multiplicity (tt) to the mean number of wounded nucleons 
(N\y)- The simple intuitive argumentation can be further quantified within 
SMES assuming the generalised Fermi-Landau initial conditions: the initial 
volume is Lorenz-contracted, V oc (v^) 1 (y/s is the c.m.s. energy of the nu- 
cleon pair), the initial energy density is given by e oc gT 4 oc (y/s — 2mi\r) ■ y/s 
(T is the initial temperature and g is an effective number of internal degrees 
of freedom at the early stage, is the nucleon mass). The pion multiplic- 
ity is proportional to the initial entropy, and the (tt)/(Nw) ratio can be thus 
calculated outside the transition region as 

(tt)/(N w ) oc V g oc g 1 ' 4 ^ - m N f 4 (JS)- l l 4 ee g 1 / 4 F . (1) 



Therefore, the (tt) / (Nw) ratio increases linearly with F outside the transition 
region, and the slope parameter is proportional to g 1 ! 4 [7]. In the transition 
region, a steepening of the pion energy dependence is predicted, because of 
activation of partonic degrees of freedom, i.e. an effective number of inter- 
nal degrees of freedom in the quark gluon plasma (QGP) is larger than in the 
hadron gas (HG): g QG p > guc- 
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The Strangeness to Pion Ratio. The energy dependence of the strangeness 
to entropy ratio is a crucial signal of the deconfinement due to its weak de- 
pendence on the assumed initial conditions. Within SMES at low collision en- 
ergies, when the confined matter is produced, the strangeness to entropy ratio 
steeply increases with collision energy. Due to a low temperature at the early 
stage (T < Tc = 170 MeV) and the high mass of the carriers of strangeness 
(ms = 500 MeV, the kaon mass) the total strangeness is oc exp(— ms/T). 
On the other hand, the total entropy is approximately oc T 3 . Therefore, the 
strangeness to pion ratio is oc exp(— ms/T) ■ T~ 3 in the HG and strongly in- 
creases with the collision energy. When the transition to a deconfined matter 
is crossed (T > Tc), the mass of the strangeness carriers is significantly re- 
duced (m s = 130 -7- 170 MeV, the strange quark mass). Due to the low mass 
(m s < T) the strangeness yield becomes (approximately) proportional to the 
entropy (both are proportional to T 3 ), and the strangeness to entropy (or pion) 
ratio becomes independent of energy in the QGP. This leads to a "jump" in the 
energy dependence from the larger value for confined matter at Tc to the value 
for deconfined matter. Thus, within the SMES, the non-monotonic energy de- 
pendence of the strangeness to entropy ratio is followed by a saturation at the 
deconfined value which is a direct consequence of the onset of deconfinement 
taking place at about 30 AGeV [1]. 

The Inverse Slope of Transverse Mass Spectra. We discuss another well 
known observable, which may be sensitive to the onset of deconfinement, the 
transverse momentum, pr, spectra of produced hadrons. It was suggested by 
Van Hove [8] more than 20 years ago to identify the deconfinement phase tran- 
sition in high energy proton-antiproton interactions with a plateau-like struc- 
ture of the average transverse momentum as a function of the hadron multiplic- 
ity 1 . According to the general concepts of the hydrodynamical approach the 
hadron multiplicity reflects the entropy, whereas the transverse hadron activity 
reflects the combined effects of temperature and collective transverse expan- 
sion. The entropy is assumed to be created at the early stage of the collision 
and is approximately constant during the hydrodynamic expansion. The multi- 
plicity is proportional to the entropy, S = s • V, where s is the entropy density 
and V is the effective volume occupied by particles. During the hydrodynamic 
expansion, s decreases and V increases with s ■ V being approximately con- 
stant. Large multiplicity at high energies means a large entropy density at 
the beginning of the expansion (and consequently a larger volume at the end). 
Large value of s at the early stage of the collisions means normally high tem- 
perature T at this stage. This, in turn, leads to an increase of the transverse 

'in the original suggestion [8] the correlation between average transverse momentum and hadron multi- 
plicity was discussed for p + p at fixed collision energy. Today we have an advantage to use central A+A 
collisions at different energies [4]. 
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hadron activity, a flattening of the transverse momentum spectrum. Therefore, 
with increasing collision energy one expects to observe an increase of both the 
hadron multiplicity and the average transverse momentum per hadron. How- 
ever, presence of the deconfinement phase transition would change this corre- 
lation. In the phase transition region, the initial entropy density (and hence the 
final hadron multiplicity) increases with the collision energy, but temperature 
T = Tc and pressure p = pc remain constant. The equation of state (EoS) 
presented in a form p(e)/e versus e shows a minimum (the 'softest point' [9]) 
at the boundary of the generalised mixed phase and the QGP. Consequently the 
shape of the pt spectrum is approximately independent of the multiplicity or 
the collision energy. The transverse expansion effect may even decrease when 
crossing the transition region [8]. Thus one expects an 'anomaly' in the energy 
dependence of the transverse hadron activity: the average transverse momen- 
tum increases with the collision energy when the early stage matter is either in 
a pure confined or in a pure deconfined phase, and it remains approximately 
constant when the matter is in a mixed phase [4, 8]. A simplified picture with 
T = Tc inside the mixed phase is changed if the created early stage matter has 
a non-zero baryonic density. It was however demonstrated [10] that the main 
qualitative features (T = const, p = const, and a minimum of the function 
p{e)je versus e) are present also in this case. In the SMES model [1], which 
correctly predicted energy dependence of the pion and strangeness yields, the 
modification of the EoS due to the deconfinement phase transition is located 
between 30 and about 160 A-GeV. Thus an anomaly in the energy dependence 
of the transverse hadron activity may be expected in this energy range. 

The energy density at the early stage increases with increasing collision en- 
ergy. At low and high energies, when a pure confined or deconfined phase is 
produced, this leads to an increase of the initial temperature and pressure. This, 
in turn, results in an increase of the transverse expansion of a matter and con- 
sequently a flattening of the transverse mass spectra of final state hadrons. The 
experimental data on the transverse mass spectra (jut = (m 2 +p\) 1 ^ 2 , m is a 
particle mass) are usually parametrized by a simple exponential dependence: 



where the inverse slope parameter T* is sensitive to both the thermal and col- 
lective motion in the transverse direction. In the parameterisation (2), the shape 
of the niT spectrum is fully determined by a single parameter, the inverse slope 
T*. In particular, the average transverse mass (%) can be expressed as: 



Hydrodynamical transverse flow with collective velocity vt modifies the Boltz- 
mann mx-spectrum of hadrons. At low transverse momenta, it leads to the 




(2) 



(m T ) = T* + m + 



(T*) 2 



(3) 



m + T* 
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result {Tki n is a kinetic freeze-out temperature): 

T&w-vr = T ^n + \mv 2 T . (4) 

A linear mass dependence (4) of T* is supported by the data for hadron spectra 
at small pp. However, for px » m the hydrodynamical transverse flow leads 
to the mass-independent blue-shifted 'temperature': 

T* - T l 1 + v T 

± high-p T ~ kln ' V I — VT ' K ' 

A simple one parameter exponential fit (2) is quite accurate up to — m = 
1 GeV for K + and K~ mesons in A+A collisions at all energies. This means 
that Tf om _ pT f» T£ igh _ pT for kaons and the energy dependence of the aver- 
age transverse mass (mr) (3) and the average transverse momentum (pr) for 
kaons is qualitatively the same as that for the parameter T*. Note that a simple 
exponential fit (2) neither works for light 7r-mesons, Tj* ow _ pT < Tj* igh _ pT , nor 
for heavy (anti)protons and (anti)lambdas, T* ow _ pT > T£ igh _ pT . This means 
that the average transverse masses, (mr), and their energy dependence for 
these hadrons are not connected to the behavior of the slope parameters in the 
simple way described by Eq. (3): one should separately consider both Tj* ow _ pT 
and Tfr igh _ pT slopes (see Refs. [11, 12] for details). 

The Dynamical Event-by-Event Fluctuations. In thermodynamics, the 
energy E, volume V and entropy S are related to each other through the EoS. 
Thus, various values of the energy of the initial equilibrium state lead to dif- 
ferent, but uniquely determined, initial entropies. When the collision energy 
is fixed, the energy, which is used to hadron production still fluctuates. These 
fluctuations of of the inelastic energy are caused by the fluctuations in the dy- 
namical process which leads to the particle production. They are called the 
dynamical energy fluctuations [5]. Clearly the dynamical energy fluctuations 
lead to the dynamical fluctuations of any macroscopic parameter of the matter, 
X, like its entropy and strangeness content. The relation between the dynami- 
cal energy fluctuations and dynamical fluctuations of the macroscopic param- 
eter X is given by the EoS. Consequently, simultaneous measurements of the 
event-by-event fluctuations of both the energy and the parameter X yield in- 
formation on the EoS. Since EoS shows an anomalous behavior in the phase 
transition region, the anomaly should be visible in the ratio of entropy to en- 
ergy fluctuations [5]. 

According to the first and the second principles of thermodynamics the en- 
tropy change SS is given as TSS = 5E+p5V. If we fix the collision geometry, 
choosing e.g. only a sample of central A+A collisions, we can expect SV = 0. 
Within SMES the ratio of entropy to energy fluctuations can be then easily 
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calculated and presented as a simple function of the p/e ratio [5]: 

Thus it is easy to predict a qualitative dependence of the R e ratio on the col- 
lision energy. Within the model, the confined matter, which is modelled as an 
ideal gas, is created at the collision early stage below the energy of 30 A-GeV. 
In this domain, the ratio p/e, and consequently the R e ratio, are approximately 
independent of the collision energy and equal about 1/3 and 0.56, respectively. 
The model assumes that the deconfinement phase-transition is of the first or- 
der. Thus, there is the mixed phase region, corresponding to the energy interval 
30-^60 A-GeV. At the end of the mixed phase the p/e ratio reaches minimum 
(the "softest point" of EoS [9]). Thus in the transition energy range the R e 
ratio increases and reaches its maximum, R e rs 0.8, at the end of the transi- 
tion domain. Further on, in the pure deconfined phase, which is represented 
by an ideal quark-gluon gas under bag pressure, the p/e ratio increases and 
approaches its asymptotic value 1/3 at the highest SPS energy 160 A-GeV. An 
estimate of entropy fluctuations can be obtained form the analysis of multiplic- 
ity fluctuations as proposed in [5]. 

At the stage of particle freeze-out, the system's entropy is related to the 
mean particle multiplicity. We assume that the multiplicity of negatively charged 
hadrons is proportional to the system entropy, S oc N^. Thus the initial en- 
tropy fluctuations are transformed into the the fluctuations of the mean mul- 
tiplicity. It is important to distinguish these dynamical fluctuations of iV_ 
formed at the initial stage of A+A reaction, from the statistical fluctuations of 
N- around iV_ at the freeze-out (see Ref. [5] for details). 

In Ref. [6] we study within the SMES the energy dependence of the dynam- 
ical strangeness fluctuations caused by the dynamical energy fluctuations. We 
define N s as a total number of strange quark-antiquark pairs created in A+A 
collision, and consider the fluctuation ratio: 

_ (8N s f/N] 

(SE) 2 /E 2 ' U) 

When T — > oo the system is in the QGP phase. Strange (anti)quarks can 
be considered as massless and the bag constant can be neglected. Then e oc 
T 4 and n s oc T 3 and consequently de/e = 4 • dT/T and dn s /n s = 3 • 
dT/T, which result in R s = (3/4) 2 0.56. In the confined phase, T < T c , 
the energy density is still approximately given by oc T 4 due to the dominant 
contributions of non-strange hadron constituents. However, the deppendence 
of the strangeness density on T is dominated in this case by the exponential 
factor, n s oc exp(— ms/T), as T << ms = 500 MeV. Therefore, at small 
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T one finds de/e = 4 • dT/T and dn s /n s = m s ■ dT/T 2 , so that the ratio 
Rs = m s/{^T) decreases with T. The strangeness density n s is small and 
goes to zero at T — > 0, but the fluctuation ratio R s (7) is large and goes to 
infinity at zero temperature limit. 

3. Signals of Deconfinement: Experimental Results 

The Pion Kink. A recent compilation of the data on the pion multiplicity in 
central Pb+Pb (Au+Au) collisions and p+p interactions is shown in Fig. 1. 

Figure 1 . The depen- 
dence of the total pion multi- 
plicity per wounded nucleon, 
(n)/{Nw), on the Fermi's 
energy measure F = (y/s — 
2m N f' A - U/?n 1/4 for cen- 
tral A+A collisions (closed 
symbols) and inelastic p + 
p(p) interactions (open sym- 
bols). The lines are calculated 
with SMES [1], The dashed 
line indicates a fit of the form 
a • F to the p + p(p) data 
(a = 1.01 ± 0.04 GeV _1/2 ). 
The solid line corresponds to 
the fit of A+A data for ini- 
tial energy above 40 A- GeV 
F ,„ „lft (a = 1.36 ± 0.03 GeV- 1/2 ). 

F (GeV ) 




One observes that the mean pion multiplicity per wounded nucleon in p + p(p) 
interactions is approximately proportional to F. For central A+A collisions, 
the dependence is more complicated and cannot be fitted by a single linear 
function. Below 40 A-GeV the ratio (tt) / (N\y) in A+A collisions is lower than 
in p + p interactions (pion suppression), while at higher energies (tt) / (N\y) is 
larger in A+A collisions than in p + p(p) interactions (pion enhancement). In 
the region between the AGS and the low SPS energy, the slope changes from 
a ^ 1.01 GeV- 1 / 2 (the dashed line in Fig. 1) to a ^ 1.36 GeV" 1 / 2 (the full 
line in Fig. 1). 

The measured increase of the slope for A+A collisions, by a factor of about 
1.3, is interpreted within the SMES as due to an increase of the effective num- 
ber of the internal degrees of freedom by a factor of (1.3) 4 = 3 and is caused 
by the creation of a transient state of deconfined matter at energies higher than 
30 A -GeV. A transition from the pion suppression [13] to pion enhancement is 
demonstrated more clearly in the insert of Fig. 1, where the difference between 
(it) I (N\y) for A+A collisions and the straight line parametrisation of the p+p 
data is plotted as a function of F up to the highest SPS energy. 
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The Strange Horn. One can argue that the strangeness to entropy ratio 
is closely proportional to the two ratios directly measured in experiments: the 
{K + )/(ir + ) ratio and the E s = ((A) + (K +K))/{tt) ratio. The energy depen- 
dence of both ratios is plotted in Fig. 2 for central Pb+Pb (Au+Au) collisions 
and p+p interactions. 
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Figure 2. The dependence of the (K + ) / {ir + ) (left) and Es (right) ratios on the collision 
energy for central A+A collisions (closed symbols) and inelastic p + p interactions (open sym- 
bols). The predictions of SMES [1] for the Es ratio are shown by a line. Different line styles 
indicate predictions in the energy domains in which the confined matter (dashed line), the mixed 
phase (dashed-dotted line) and the deconfined matter (dotted line) are created at the early stage 
of the A+A collision. 

For p+p interactions both ratios show monotonic increase with energy. How- 
ever, very different behavior is observed for central Pb+Pb (Au+Au) collisions. 
The steep threshold rise of the ratio characteristic for confined matter then set- 
tles into saturations at the level expected for deconfined matter. In the transition 
region (at low SPS energies) a sharp maximum is observed caused by a higher 
strangeness to entropy ratio in the confiend matter than in the deconfiend mat- 
ter. As seen in Fig. 2 the measured dependence is consistent with that expected 
within the SMES. 

The Step in Slopes. The energy dependence of the inverse slope parameter 
fitted to the K + and K~ transverse mass spectra for central Pb+Pb (Au+Au) 
collisions is shown in Fig. 3 [4]. The striking features of the data can be sum- 
marised and interpreted as follows. The T* parameter increases strongly with 
collision energy up to the lowest (30 A-GeV) SPS energy point. This is an 
energy region where the creation of confined matter at the early stage of the 
collisions is expected. Increasing collision energy leads to an increase of the 
early stage temperature and pressure. Consequently the transverse activity of 
produced hadrons, measured by the inverse slope parameter, increases with 
increasing energy. The T* parameter is approximately independent of the col- 
lision energy in the SPS energy range. In this energy region the transition 
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between confined and deconfined matter is expected to be located. The re- 
sulting modification of the equation of state "suppresses" the hydrodynamical 
transverse expansion and leads to the observed plateau structure in the energy 
dependence of the T* parameter. At higher energies (RHIC data), T* again 
increases with the collision energy. The equation of state at the early stage 
becomes again stiff, the early stage temperature and pressure increase with 
collision energy, and this results in increase of T* too. 
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Figure 3. The energy dependence of the inverse slope parameter T* for K + (left) and K~ 
(right) mesons (see Ref. [4]) produced at mid-rapidity in central Pb+Pb (Au+Au) collisions at 
AGS (triangles), SPS (squares) and RHIC (circles) energies. 



Kaons are the best and unique particles among measured hadron species for 
observing the effect of the modification of the equation of state due to the onset 
of the deconfinement in hadron transverse momentum spectra. The arguments 
are the following. 

1) . The kaons m-p-spectra are only weakly affected by the hadron re- 
scattering and resonance decays during the post-hydrodynamic hadron cascade 
at the SPS and RHIC energies [11]. 

2) . A simple one parameter exponential fit (3) is quite accurate for kaons in 
central A+A collisions at all energies. This simplifies strongly an analysis of 
the experimental data. 

3) . The high quality data on mr-spectra of K + and K~ mesons in central 
Pb+Pb (Au+Au) collisions are available in the full range of relevant energies. 

The Dynamical Event-by-Event Fluctuations. The ratios of entropy to en- 
ergy (6) and strangeness to energy (7) dynamical fluctuations calculated within 
SMES are presented in Fig. 4. We find a non-monotonic energy dependence 
of R e with a maximum at the boundary between the mixed phase and the QGP 
[5]. A pronounced minimum-structure is expected in the dependence of R s on 
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the collision energy [6]. It is located at 30 60 AGeV, where the mixed phase 
is created at the early stage of A+A collision. 
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Figure 4. The collision energy dependence of the relative entropy to energy fluctuations 
R e (6), left (see Ref. [5]), and strangeness to energy fluctuations R„ (7), right (see Ref. [6]), 
calculated within SMES. 



The experimental data on the energy dependnece of R e and R s ratio are not 
yet available. Both entropy and strangeness fluctuation measures, R e and R s , 
show anomalous behavior in the transition region: the maximum is expected 
for R e and the minimum for R s . Consequetly, even a stronger anomaly is 
predicted for the ratio: 

' Re (*JV_) 2 /Fi 

Experimental measurements of R s / e may be easier than the measurements of 
R e and R s because the ratio R s i e requires measurements of particle multi- 
plicities only, whereas both R e and R s involve also measurements of particle 
energies. 

4. Conclusions 

The energy scan program at the CERN SPS together with the measure- 
ments at lower (LBL, JINR, SIS, BNL AGS) and higher (BNL RHIC) energies 
yielded systematic data on energy dependence of hadron production in cen- 
tral Pb+Pb (Au+Au) collisions. Predicted signals of the deconnnement phase 
transition, namely anomalies in the energy dependence of hadron production - 
the pion kink [1, 7], strange horn [1] and the step in slope of i\T-mesons [4] - 
are simultaneously observed in the same domain of the low SPS energies and 
presented in the left panel of Fig. 5. 
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Figure 5 . The collision 
energy dependence of 
{n)/{N w ),{K+)/(-K+),md 
T*(K + ). The predicted sig- 
nature of the deocnfinement 
transition - the pion kink 
[1,7], the strange horn [1] and 
the step in mr-slopes [4] are 
observed simultaneously in 
the same domain of collision 
energies F = 2 -f- 3 GeV 1/2 . 
The open symbols correspond 
to the data in p + p collisions. 
The right panel shows the 
prediction of SMES for 
R s /e (8) constructed from 
fluctuations of strange and 
non-strange hadron multi- 
plicities [6]. A pronounced 
minimum of R s / e is predicted 
in the collision energy domain 
in which deconfinement 
transition is located. 
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The anomalies in the energy dependence of the hadron production are seen in 
central A+A collisions and absent in the data of p + p reactions. They indi- 
cate that the onset of the deconfinement in central Pb+Pb collisions is located 
at about 30 A-GeV. The theoretical picture is however far from being com- 
plete. It seems that the Landau-type initial conditions are not appropriated to 
describe both the multiplicity and transverse momentum data at the same time. 
The initial conditions with spatially spread energy and longitudinal velocity 
distributions are needed to reproduce all main features of the data simultane- 
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ously [14]. The analysis of the data from the energy scan program is also still 
in progress. In particular, first results at 20 A-GeV are soon expected. We hope 
that the properly analysed event-by-event fluctuations (see Ref. [5]) may also 
be sensitive to the onset of the deconfinement. Especially promising looks a 
new measure R s / e (8) constructed from the fluctuations of strange and non- 
strange hadron multiplicities [6] and shown in the right panel of Fig. 5. 
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